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ABSTRACT 

We have obtained spectra with resolution ^2.5 A in the region w 7500-9500 A for 116 
red giants in 5 Galactic globular clusters and 6 old open clusters (5 with published 
metallicities, and 1 previously unmeasured). The signal-to-noise ranges from 20 < S/N 
< 85. We measure the equivalent widths of the infrared Ca II triplet absorption lines in 
each star and compare to cluster metallicities taken from the literature. With globular 
cluster abundances on the Carretta & Gratton scale, and open cluster abundances 
taken from the compilation of Friel and collaborators, we find a linear relation between 
[Fe/H] and Ca II line strength spanning the range —2 < [Fe/H] < —0.2 and 2.5 < 
(age/Gyr) < 13. The reference abundance scales appear to be consistent with each 
other at the ~0.1 dex level. Alternate choices for metallicity scales can introduce 
curvature into the relation between [Fe/H] and Ca II equivalent width. No evidence 
for an age effect on the metallicity calibration is observed. Using this calibration, we 
find the metallicity of the massive, old, open cluster Trumpler 5 to be [Fe/H] = —0.56 
±0.11. This is the first spectroscopic abundance measurement for Trumpler 5, and is 
lower by ps 0.3 dex than estimates based on the cluster colour-magnitude diagram. 
Considering the 10 clusters of known metallicity shifted to a common distance and 
reddening, we find that the additional error introduced by the variation of horizontal 
branch/red clump magnitude with metallicity and age is of order ±0.05 dex, which 
can be neglected in comparison to the intrinsic scatter in the method. The results are 
discussed in the context of abundance determinations for red giants in Local Group 
galaxies. 

Key words: techniques: spectroscopic - stars: abundances - stars: late-type - open 
clusters: individual: Trumpler 5 



1 INTRODUCTION 

One of the most widely applied techniques for the deriva- 
tion of abundances of individual red giants in very dis- 
tant systems (globular clusters or external galaxies) is the 
measurement of the strength of the near-infrared lines of 
Ca II at AA = 8498.02 8542.09, 8662.14 A . As noted 
in pioneering papers by lArni androff fc ZinrJ (ll988T> and 
lArmandroff fc Da Costal il99 it) , the calcium triplet (CaT) 
is attractive for several reasons: the red giants are near their 
brightest in this wavelength range; the lines are so broad 
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that they can be measured at moderate resolution; they 
arise between excited states of the Ca + ion and so are rel- 
atively unaffected by interstellar absorption; and they lie 
in a window of small telluric absorption between H2O and 
O2 bands. Ca II is the dominant ionization stage in red gi- 
ant branch (RGB) stars, and because the continuous opac- 
ity in the photospheres is dominated by the H~ ion, the 
relation between line s trength and surface g ravity is theo- 
retically well-modelled ll0rgensen et al.ll992f) . Additionally, 
the CaT is attractive because calcium abundances are not 
expected to be affected by proton-capture nucleosynthesis 
in the intermediate- and low-mass stars of the RGB. This 
distinguishes it from lighter elem ents such as O, Na, Mg, 
Al, and Si (e.g.. llvans et al]l200lT) . Thus the amount of cal- 
cium in an RGB envelope should be representative of the 
primordial abundances in the star. 
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However, the matter is complicated by the very strength 
of the lines that makes them attractive in the first place. 
The CaT lines are strongly saturated and extremely broad. 
The line strength thus depends strongly on the tempera- 
ture structure of the upper photosphere and chromosphere, 
regions which are notoriously difficult to model in red gi- 
ants. Ca is an a-element, so its abundance relative to the 
iron-peak elements is expected to vary with environment 
based on the relative importance of Type I and Type II su- 
pernovae to chemical evolution. Quantifying the impact of 
[a/Fe] variation on abundances derived from the CaT is non- 
trivial, because the a-enhancement affects both the stellar 
temperature and the electron density in the atmosphere. A 
detailed theoretical understanding of the relation between 
CaT line strength, metallicity, and the evolutionary state of 
a giant star is quite challenging. 

In contrast, a large body of empirical work shows that 
the CaT line stren gth is highly sensitiv e to overall stel- 
lar metallicity (e.g., Rutledee et aljlf 997bl hereafter R97b). 
Their calibration uses a simple linear relation between the 
CaT pseudo-equivalent widths and V magnitudes of indi- 
vidual red giants to effectively remove the T e ff and log g 
dependencies. A spectrophotometric index that is sensitive 
only to metallicity is thus defined and calibrated to a sam- 
ple of known metallicity. Qualitative theoretical understand- 
ing of this line strength behavior is based on the response 
to changes induced in the stellar atmosphe re structure by 
metallicity variations llj0rgensen et alJ ll992l, 

Because the CaT index is empirically calibrated, it is 
potentially dangerous to apply the technique to stellar pop- 
ulations beyond the range of the calibration sample. The 
m ost complete such calibr ation to date has been published 
bv lRutledge et alJ (|l997al . hereafter R97a) and R97b, who 
found a linear relation between the reduced CaT equiva- 
lent widths and metallicities o n the high-dispersion glob- 
ular cluster abundance scale of ICarretta fc Crattonl l)l997i 
hereafter CG97). This calibration, while extremely useful, 
only covers the metallicity range from —2.1 < [Fe/H] < 
—0.6 and the age range from 9 < (age/Gyr) < 13. This is 
currently one of the limiting uncertainties in understand- 
ing the age- metallicity relation and ab undance gradient 
of the Large Magellanic Cloud (LMC) jCole et alJ l2000t 
ISmecker-Hane et alJ 2003). for example, and might be prob- 
lematic for planned studies of mixed-age and metal-rich sys- 
tems like the M3I halo. 

It is important to develop a detailed understanding of 
the calcium triplet technique because at present it is the 
only way to efficiently make spectroscopic abundance as- 
sessments of red giants in galaxies beyond the Milky Way 
halo; this makes it an essential observational tool for the 
understanding of chemical evolution of galaxies. Here, we 
show that the age effect is a negligible source of error for 
metallicities derived from the CaT index. One of us (TLB) 
will expand on this result in a forthcoming analysis of high- 
and low-dispersion spectra of numerous Galactic star cluster 
giants, producing a better understanding of the relationship 
between the calcium abundance and the CaT line strength 
as part of her Ph. D. thesis. 

In this paper we move towards the calibration of the 
CaT index as a metallicity indicator for individual red gi- 
ants at younger ages than the globular clusters that de- 
fine the present CaT abundance scale. The following section 



describes the sample of observed clusters, the data reduc- 
tion process, and the measurement of radial velocities and 
equivalent widths. Section discusses the relation between 
calcium line strength and metallicity for the 10 clusters of 
known abundance, and the effects of age on the metallicity 
measurement. In §Z\ we apply our calibration of the CaT 
index to the little-studied massive cluster Trumpler 5, pro- 
ducing the first spectroscopic metallicity assessment for the 
cluster. We consider the systematic errors introduced by tak- 
ing a stellar sample of mixed age and metallicity, as in Local 
Group galaxies in section |H| Lastly, we briefly discuss the 
choice of reference metallicity scales and their effects on our 
calibration. 



2 DATA: OBSERVATIONS & REDUCTIONS 

The data were obtained under ESO programme 70.B-0398, 
a study of the metallicity distribution function of the cen- 
tral bar of the LMC. Because most LMC field stars are 
much younger than the Galacti c globular clusters feutcherl 
If 977tlSmecker-Hane et alJl2002T) and have metal licities near 
the upper end of the R9 7b calibration ijCole et al.l 2000; 
ISmecker-Hane et alJ EooS) . we observed a sample of five 
globular clusters and six old open clusters to calibrate our 
LMC sample and to begin the extension of the CaT abun- 
dance scale to younger ages and higher metallicities. The 
observed clusters and their relevant physical parameters are 
summarized in Table Q 



2.1 Target Selection 



We 



chose to 
been 



observe 



globular 
used 



clusters 
in the 



that have al- 
calibration of 



the CaT (e.g. 


Smecker-Hane et al.l 120031 iTolstov et al 


2001 
1995 




200ft iR.utleoVe et alJll997al: iGeisler et al 


: iDa Costa & Armandroff 1995: Suntzcff et al. 


1993: 


Armandroff & Da Costal If 99ll lOlszewski et alJlf99I 


), lim- 



ited to those with good late-December visibility from north- 
ern Chile. All of these clusters, except NGC f85f, have 
had high-dispersion metallicity analyses published by CG97; 
NGC f85f has been measuered with the CaT technique 
and calibrated to the CG97 scale by R97b. While the abso- 
lute ages of globular clusters remain uncer tain, our sample 
has be en placed onto a consistent scale bv lSalaris fc Weissl 
(2002). Table |2| shows the individual red giants observed 
in each cluster, with the adopted V magnitude, derived ra- 
dial velocities and CaT equivalent widths, and references to 
the identifications and photometry. The number of stars ob- 
served in each cluster depends on the cluster angular size 
and density and the spectrograph field geometry, hence rel- 
atively few stars have been observed in nearby (47 Tuc) or 
sparse (NGC 2298) clusters. 

Open cluster targets were chosen from the iFriel et al] 
(2002, hereafter F02) list of old clusters with metallicity es- 
timates based on medium-resolution measurements of iron 
and iron-peak-element line strengths. F02 also include age 
estimates for each cluster in their sample based on the 
colour-magnitude diagram morphology of the clusters. We 
chose to adopt this relatively homogeneous set of measure- 
ments as our reference sample for open clusters. The true 
values of age and metallicity for many clusters are surely still 
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systematically uncertain by substantial amounts; we adopt 
the numbers from Friel and collaborators here because they 
place a large number of clusters onto a consistent scale that 
is compatible with the CG97 globular cluster scale (see be- 
low). 

We picked specific equatorial and southern clusters to 
observe based on richness and coverage of the age-metallicity 
plane. NGC 2682 (M67) is a much-observed target, although 
it is not a convenient one for FORS2 because of its extremely 
bright member stars and large angular size. However, it re- 
mains a useful calibrator for metallicities close to solar. Our 
open cluster RGB targets are listed in Table El with refer- 
ences to the photometry and identifications. 

2.1.1 Star Cluster Metallicity Scales 

There is significant uncertainty in the absolute metallicity 
scales of globular and open clusters. Efforts to establish large 
samples of clusters analyzed in a uniform way (consistent 
effective temperature scales, model atmospheres, and lines 
analyzed) have not generally overlapped between the two 
cluster families, ft is therefore important to establish that 
our adopted cluster metallicities are not grossly incompati- 
ble with each other. 

The F02 spectroscopic indices were calibrated to a sam- 
ple of field and cluster stars with high-dispersion analyses. 
Two globular clusters, M71 and M5, were incl uded; the ref- 
erenc e metalliciti es were taken to b e —0.79 dSneden et alJ 
Il994l) and -1.17 JSneden et al.ll992l) respectively. These are 
lower by 2a than the metallicities on the CG97 scale, —0.73 
±0.03 and -1.12 ± 0.03 (R97b). Mos t of F02's calibra- 
tion sa mple were field giant s measured bylCottreh^J^jjneder 
( 1986) and analyzed by IShetrone. Sneden fe Pilachowski 
( 1993) . These authors used MARCS atmosphere models 
iGustafsson et al.lll975l) and effective temperatures derived 
from imposing excitation equilibrium of Fe I and Ti I in 
the spectra. I n con trast, CG97 used stellar atmospheres 
from iKuruczl 1^993) and derived T e / / from a V— K color- 
temperature relation. 

While the two halves of our reference sample are 
inhomogeneous, there is good reason to think they are 
not inconsistent. Inspection of abundances derived from 
the same observational material using MARCS versus Ku- 
rucz atmospheres (e.g. iKraft fe Ivansl 12003) indicates that 
the MARCS-based abundances tend to be lower by ~ 
0.08 dex. Comparing the differences due to T e ff scales 
is more problematic, given the extremely heterogeneous 
nature of the material. We can estimate the effects of 
the inconsistent T e // scales by noting that: a) temper- 
atures derived from the infrared flux method (IRFM; 
Alon so. Arriba s & Mar tinez- Rogeilll999D tend to be in ex- 
cellent agreement with those deri ved from Fe I ex citation 
equilibria in MARCS models (e.g.. lHill et al.ll200ol) . and b) 
IKraft fc Ivarisl (12003ft find that IRFM temperatures tend to 
be ~50 K lower than color temperatures, resulting in abun- 
dance values higher by « 0.05 dex when the former scale is 
used. 

There are two implications of these trends: first, that 
the systematic offsets in abundance between the adopted 
globular and open cluster metallicity scales are comparable 
to or smaller than the internal uncertainties within the re- 
spective scales; and second, that when comparing the F02 



metallicities to the CG97 metallicities, the offset caused by 
choice of model atmospheres may be of similar magnitude 
but opposite sign to the offset caused by choice of T e ff 
scales. However, offsets between the two metallicity scales 
of order 0.1 dex are probable; the effects of adopting differ- 
ent scales are discussed further in section |H] 



2.1.2 Trumpler 5 

In addition to the five well-studied open clusters, we ob- 
served the important, but understudied, cluster Trumpler 5 
(= Collinder 105 = C 0634+094) in order to provide a test 
of our abundance calibration. Tr 5 is among the most mas- 
sive open clusters in the Galaxy (Kaluzny 1998j, hereafter 
K98), but its high reddening has made it an unattractive 
object of study for the most part. Tr 5 is an interesting ob- 
ject because there is no published spectroscopic abundance 
measurement, so we can immediately apply our calibration 
to a derivation of its metallicity. One photometric abun- 
dance estimate places it near solar metallicity (K98), and 
it is of a richness and distance that would make it a more 
convenient future calibrator of the CaT in the moderate- 
abundance regime than M67 has been. It will be critical 
to extend this calibration to solar abundance and above 
as the CaT method is applied throughout the giant and 
dwarf galaxies of the Local Group. The metallicity of Tr 5 
is in itself a valuable datum, because of the possible mem- 
bers hip of the N-type carbon star V493 Mon in the clus- 
ter llKalinowski. Burkhead. fc Honevcutd Il974l) . An accu- 
rate measure of the cluster metallicity would in turn lead 
to a more precise age estimate for the cluster, and possi- 
bly to new constraints on models of carbon star formation. 
We selected can didate members of Tr 5 using 2MASS J and 
Ks photometry JSkrutskie et al.lll997ft and later made the 
cross-identification with BV photometry of K98; targets are 
identified in Table 2] Note that three stars appeared to be 
differentially reddened in the (B— V, V) CMD, and that lack- 
ing prior radial velocity information to guide selection, two 
foreground stars were included in the sample. V493 Mon was 
observed in order to compare its radial velocity to that of 
the cluster (no CaT abundance estimate can be made for 
C stars). 

2.2 Data Acquisition and Analysis 

The observations were made in Visitor Mode at the Yepun 
(VLT-UT4) 8.2-metre telescope of ESO's Paranal Ob- 
servatory, with the FORS2 spectrograph in multiobject 
(MOS) mode, between 24 and 26 December 2002. We used 
the 1028z+29 volume-phased holographic grism with the 
OG590+32 order blocking filter. In MOS mode, the FORS2 
field is covered by an assembly of 19 mechanical slit jaws, 
20"-22" long, that can each be arbitrarily positioned along 
the horizontal axis of the field (East-West for instrument 
rotator angle = 0°). Throughout, we used a slit width of 
1". The red-optimized MIT/LL 2kx2k pair of CCDs were 
in use, with readnoise of 2.7 electrons and gain of 0.8 e~ 
ADU -1 . In the standard configuration, the CCD pixels are 
binned 2x2, giving a plate scale of 0.25" per pixel in the 
spatial direction. This setup yields spectra covering « 1700 
A, a central wavelength near 8500 A and a dispersion of 
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« 0.85 A pix -1 (resolution 2-3 A). The FORS2 field is 6.'8 
square, but it is limited to 4'.8 useful width in the dispersion 
direction in order to keep the wavelength range of interest 
from falling off the edges of the CCD. The typical exposure 
time per cluster was 3x40 seconds, although the presence of 
9th magnitude stars in M67 led us to use exposures of just 
3 seconds. Skies were clear during the entire run; the seeing 
during these observations varied between O'.'S < FWHM < 
1'.'2, with a median stellar FWHM = 0'.'6. 

Astrometry of the targets was obtained via short pre- 
imaging exposures with FORS2, taken several weeks before 
the observing run. The relative positions of the slit jaws were 
precisely set before the observing run using the FORS In- 
strument Mask Simulator (FIMS) software tool distributed 
by ESO. Through-slit images were then taken immediately 
prior to the science exposures to confirm accurate slit cen- 
treing, and the telescope was offset by up to ~ 0.25" (1/4 
of the slit width) when necessary. In cases where the targets 
spanned a large range of x- values on the chip, it was noted 
at the time of observation that stars at one end of the range 
could not be perfectly centred simultaneously with stars on 
the opposite end. The differential offset was later estimated 
to amount to roughly half a pixel over a 750 pixel baseline. 
Implications of this are briefly discussed in the following sec- 
tion. 

Calibration exposures were taken during the daytime 
under the FORS2 Instrument Team's standard calibration 
plan. These consisted of lamp flatfields with two different il- 
lumination settings and a helium-neon-argon arc lamp expo- 
sure. Two lamp settings are required for FORS2 flatfielding 
to correct for extra reflections introduced into the optical 
path of the calibration assembly (T. Szeifert 2003, private 
communication). Because of the large number of setups in 
our programme, twilight sky flats were deemed impractical. 
The basic data reduction was performed in IRAF 1 . We fit 
and subtracted the scaled overscan region, trimmed the im- 
age, and divided by the appropriately combined lamp flats 
within IRAF's ccdred package. 

The spectroscopic extractions were done using hydra, 
an IRAF package for handling multislit spectra. Our tar- 
gets were bright enough that the trace could be extracted 
directly from the science exposures. Across the field of view, 
the curvature of the spectral trace varied significantly, but 
could in all cases be well-fit with a low-order polynomial. 
Because of the potential for spectrograph flexure during 
the night and the high probability of small, uncorrected 
slit-c entreing errors, we used n ight sky emission lines of 
OH JOsterbrock fe Marte]Hl992T) and 2 JOsterbrock et alJ 
to dispersion-correct the spectra instead of the day- 
time arc lamp exposures. Typically, 25-35 strong lines were 
used in the wavelength solutions, which had characteristic 
root-mean-square (r.m.s.) scatter of 0.04-0.08 A. The ex- 
ception to this rule was M67, for which the exposures were 
so short that no sky lines were visible. The spectral win- 
dow covered depends on the position of the slit along the 
dispersion direction of the CCD; some spectra reach as far 
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into the blue as 7300 A, and some as far into the infrared as 
10100 A, but most are centred near 8500 A and cover the 
approximate range 7600 A < A < 9350 A. 

Extraction to I-dimensional spectra was done within 
the apall tasks. Sky subtraction was performed using 1- 
dimensional fits across the dispersion direction. Because the 
targets are bright compared to the sky and the slitlets are 
long compared to the seeing disc, this presented few difficul- 
ties, except when the stars fell near the ends of the slitlets. 
In these cases, the sky region was chosen interactively and 
adjusted to produce the cleanest extracted spectrum in the 
region around the CaT. Dispersion-corrected spectra were 
combined using scombine to reduce noise from cosmic rays 
and CCD defects. The stars were continuum normalized by 
fitting a polynomial to the spectrum, excluding the Ca II 
lines and regions of strong telluric absorption. We achieved 
typical signal-to-nosie ratios of 20 < S/N < 85 per pixel. 
Sample spectra showing the region of the CaT appear in 
Figure 

2.3 Radial Velocities 

Stellar radial velocities are important for deriving the ex- 
pected centroid of the Ca lines in the equivalent width mea- 
surement, and in some cases for establishing cluster member- 
ship. We are not primarily interested in the absolute radial 
velocities of the target stars, so we did not observe radial 
velocity standards. Ra ther, we cross-correlated our spectra 
jTonrv fc Davis! I1979T) with template stars already in our 
abundance programme. This guarantees a good match in 
spectral type between the templates and the stars with un- 
known velocities. The templates were picked based on re- 
liable velocity measurements from the literature. We chose 
measureme nts of 12 stars in NGC 4590, NGC 2298, and 
NGC 1904 dCeisler et alJll995l) . 10 stars in Berkeley 20 and 
Berkeley 39 jFriel et all 120021) . and 2 stars in Melotte 66 
jFriel fc Janeslll993l) . We used the IRAF f xcor task to per- 
form the cross-correlation for each of the 24 templates, and 
determined the stellar radial velocities using the average of 
template offsets weighted by the error in the correlation re- 
sult and the height of the correlation peak. The velocities 
were corrected to the heliocentric reference frame within 
fxcor, after suitable alteration of the ESO image headers 
to the format demanded by the task. 

When the stellar image is smaller than the slit width, 
the random error in the velocities can be much smaller than 
the systematic error due to slight misalig nments between 
the s lit centre and the stellar centroid fe.g.. llrwin fc Tolstovl 
2002). When the seeing was better than 1", we evaluated the 
slit-centreing accuracy of each target by visual inspection of 
the through-slit images. Where necessary, we applied cor- 
rections to the velocities of both templates and targets. Our 
spectral resolution corresponds to ~ 29.5 km sec -1 pix -1 , 
yielding corrections that range from |Av| = 0-32 km sec -1 . 
In most cases, the offset between slit and star was strongly 
correlated with the stellar a>position on the CCD. We es- 
timate that the limit of centroiding accuracy is roughly a 
quarter of a pixel, or w 7 km sec -1 . We therefore adopt this 
value as a floor in the error estimates on the heliocentric ra- 
dial velocities in Tables|21Bl This is a conservative estimate, 
as the computed velocity dispersions for some of the clusters 
are less than this floor. 
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The mean derived cluster velocities are given in Table 
with the r.m.s. dispersion and a previous measurement 
from the literature, where available. In most cases, the clus- 
ter velocities are in reasonable agreement with previously 
published values. For the globular clusters, there is a mean 
velocity shift of 4.7 km sec -1 compared to the literature ve- 
locities. Exposure times for M67 were so short that no sky 
lines could be measured for a good dispersion solution, so we 
simply adopted the shape of the function from other targets, 
and shifted the zeropoint using the telluric A band of O2. 
Therefore we have no independe nt measure of V r for M 67; 
radial velocities from the list of iMathieu et all lll986fi are 
listed in Table 01 For comparison, the mean V r from F02 is 
shown in Table [K] 

No measurement of Trumpler 5's radial velocity could 
be found in the literature; it will be discussed further in 
section [I] Of the four remaining open clusters, agreement 
with previously published radial velocities is relatively good 
for three: Be 20, Be 39, and Mel 66 are measured within 
10 km sec -1 of the published values. However, NGC 2141 is 
found to be blueshi fted more than 30 km sec -1 from the ve- 
locity reported in lFriel. Liu, fc Janea <ll989f) . The source of 
the discrepancy is unknown, but our value is in good agree- 
ment with recent medium-dispersion measurements inde- 
pendently made at the Lick Observatory by one of us (TLB). 
In all, five stars with velocities that were highly discrepant 
from their respective cluster means were judged nonmem- 
bers and excluded from further analysis; they are flagged in 
the appropriate tables. 



2.4 Equivalent Widths 

We used interactive software kindly provided by G. Da Costa 
(1999, private communication to TSH), and later modi- 
fied by A AC, to fit the line strengths of the Ca II lines. 
The line and cont i nuum bandpasses are as defined in 
lArmandroff fc Zinnl jl98STl . Doppler-shifted in each case to 
the measured radial velocity of the target. The CaT lines 
can be contaminated by weak neutral metal lines, and the 
continuum may be affected by weak molecular bands. Both 
effects make it impossible to measure the "true" equiv- 
alent widths of the CaT lines in isolation. The spectra 
need not be flux-calibrated because the continuum slopes 
of red giants in this wavelength range are generally close 
to flat, so a simple normalization by an arbitrary polyno- 
mial suffices f or measurement of the so-c alled "pseudoequiv- 
alent width" jArmandroff fc Zinnlll98l) . The measured line 
strengths depend on spectral resolution, and measurements 
of a given star by different groups ofte n dev iate from each 
other. However, iRutledee et al.l l^)97a) and ICenarro et alJ 
( 2001aJ) have shown that measurements of the CaT index on 
various systems are generally readily transformable, so this 
is not a worry. 



2.4-1 Fitting the Line Profiles 

Our initia l approach, followin g lArmandroff fc Da Costal 
Jl99ll) and lSuntzeff et alJ lll993l) , was to fit a Gaussian func- 
tion to the Ca II lines. However, it soon became apparent 
that for the stars in the high-metallicity half of our sam- 
ple, the line profiles strongly deviate from Gaussian because 



of the strong damping wings. This effect is even stronger 
for stars in our LMC sample (Cole et al. 2003, in prepa- 
ration), because of their lower average surface gravities and 
hence broader lines. Therefore we modified the profile-fitting 
software to fit each line with the sum of a Gaussian and a 
Lorentzian function: 



F{\) = C(A) 
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(1) 



where the flux at wavelength A is equal to the continuum 
level C(A) minus the summed Gaussian and Lorentzian func- 
tions, with respective amplitudes ct\ and 0:2, and common 
line centre Ao. The width parameters a and F have their 
usual definitions. Note that the continuum level is effectively 
fixed at 1 by our normalization procedure. 

The best-fitting parameters were determined in a least- 
squares way, usi ng the Leven b erg-M arquardt algorithm as 
implemented by IPress et alJ dl992h . Initial guesses were 
based on the line depth and breadth in order to maximize 
speed of convergence. There are no significant differences be- 
tween the shape of the composite profile and the actual line 
shapes, a much more satisfactory situation than the Gaus- 
sian fits. The fits are then integrated over the line bandpasses 
to give the equivalent widths. Error estimates are obtained 
by measuring the root-mean-square scatter of the data about 
the fit function. 

The composite function fit provides a more satisfactory 
fit to the line shape, but does it change the resulting equiv- 
alent width measurement? We tested the effect of choice of 
fitting function by remeasuring our targets with the ampli- 
tude of the Lorentzian component equal to zero. Figure |21 
shows the comparison. The relation between the Gaussian- 
fit equivalent width, EWg, and the composite-fit, EWg+l 
is linear over much of its range, but there is a change in 
slope at EWg+i ~ 6.5 A. A least-squares fit to the globular 
cluster stars produces the relation 



YW G = 0.940 EWg+l ~ 0.241, 



(2) 



with r.m.s. scatter of 0.11 A about the fit for the 57 globular 
cluster stars. This relation is plotted with a dashed line in 
the upper panel of Fig. [5] The bottom panel of the figure 
shows the residuals to this fit, demonstrating that the open 
cluster stars have similar scatter to the globulars, but are 
systematically offset by a mean residual of —0.2 A (~ 2- 
3%) from the fit. This can be intuitively understood by not- 
ing that with increasing line strength, an increasingly larger 
fraction of the absorption occurs in the line wings, resulting 
in an ever-larger fraction of absorption missed by the Gaus- 
sian fit. The offset corresponds to a difference of 0.05-0.15 
dex in the derived metallicities of strong-lined stars. 

The "gross underestimation" of the strength of the line 
wings by a Gaussian fit was already noted by R97a. Because 
of the increasingly important contribution of the damping 
wings to the total line width for the strongest lines, we there- 
fore consider only the equivalent widths derived from com- 
posite line fits in the remainder of this paper. In the next 
section, we show that with this choice, the CaT index main- 
tains sensitivity over the range of metallicities considered 
here. However, the composite fit is noisier than the Gaus- 
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sian fit because the line depth is small in the far wings. By 
adding artificial noise to one of our spectra and remeasuring 
the equivalent widths, we estimate that for signal-to-nosie 
ratios less than » 15-20, there is no significant difference 
between the two profile shapes. Thus where spectra of low 
S/N and/or of low metallicity stars are concerned, a fit with 
fewer free parameters, such as a Moffat function (R97a) may 
be preferable. 

2.4-2 Forming the CaT Index 

Metallicity measurements using the CaT index use a lin- 
ear combination of the thre e individual Ca II lines to form 
the line strength index EW dArmandroff fc Da Costall99lT) . 
Various authors have used different procedures to form 
this index, depending on spectral resolution and signal- 
to-noise. Some aut hors have excluded the weakest (A = 
8498 A) line (e.g.. ISuntzeff et al.lll993 ICole et al]l200Cft . 
Others have used all thr ee lines, either weighted (R97a, 
[gme cker-Hane et alJl20o"i?) or unweighted (Ols zewski et alj 
1991). Because of our extremely high signal-to-noise ratios, 
we use an unweighted sum of the three lines. The fractional 
errors on the individal lines are roughly equal, and variations 
within the sample make it inadvisable to apply a uniform 
weighting scheme, uniform weig hting scheme for all clusters . 
It has been shown by R97b and ISmecker-Hane et al.l (I2003T) 
that the various means of creating EW produce well-defined 
linear mappings between the resulting index values. For the 
remainder of this paper, we adopt the convention 

EW = WW + V^8542 + W m &2 ■ (3) 

Stellar effective temperature and surface gravity play a 
strong role in forming the Ca II lines in red giants, and so 
the metallicity dependence is typi cally found using the so- 
called "reduced" equivalent width iArmandroff fc Da Costal 
Il99ll) . This quantity is defined as 

W' = EW±/3(V- Vhb), (4) 

where the definition exploits the empirical fact that in an 
individual globular cluster, all red giants at a given V mag- 
nitude have the same EW, and that EW (V) increases lin- 
early up the RGB. Vhb here is the mean magnitude of the 
horizontal branch; its introduction removes any dependence 
on the cluster distance or redd ening. Other paramet eriza- 
tions for W' are possible (e.g.. lolszewski et al.lll99ll) . but 
Equation 0] has the advantages that a huge amount of V 
band photometry is available in the literature, and that no 
additional color-dependence is introduced. R97a calibrated 
Equation fusing data for 52 Galactic globular clusters and 
derived a slope j3 = 0.64 ±0.02 A mag , independently 
of metallicity, across the range —2.1 < [Fe/H] < —0.6. 
Consistency with this constant /3 has been found in smaller 
studies incl uding more metal-ric h (open) clu sters, such as 
Melo tte 66 iOlszewski et alJl99ll). NGC 2477 (iGeisler et al] 
[ 1995D. and M67 JOlszewski et al.l Il99il ICole et all l2000t 
ISmecker-Hane et alJ 12003^ . Defined as in Equation 21 ev- 
ery cluster can be assigned a single value of W' that is 
strongly correlated with metallicity. The relation is linear for 
the CG97 metallicity scale, but mildly nonlinear for other 
choices (R97b). 

Because we use a different functional form to fit the 
Ca II lines of our stars, we also need to rederive f3. We did so 



by deriving the slope for each cluster and taking the average, 
weighted by the error in the individual values. The resulting 
slope is p = 0.73 ±0.04. This value is higher than found for 
the globular cluster-only sample of R97a. For reference, the 
slope determined from our sample of five globulars alone 
is Pgiob = 0.66 ±0.03, which is indistinguishable from the 
favored value of R97a, but <2<r di fferent from our value fo r 
all 10 clusters. From the models of IJeirgensen et al] lll992l) . 
it can be inferred that an increase in j3 with metallicity is 
in qualitative agreement with theory; direct comparison is 
difficult because of the interplay between abundance, gravity 
and temperature and how these quantities vary with age 
among RGB stars. 

The equivalent width measurements for the individual 
stars are shown in Figure|3] with photometry drawn from the 
sources listed in Tables|5]&|3] Error bars are omitted so that 
the best-fitting line (with slope = 0.73) to each cluster can 
be clearly seen. Note that there is significant scatter about 
the cluster fiducials even for the well-studied globulars. The 
mean r.m.s. scatter about each fiducial is ±0.24 A, larger 
than our estimated measurement errors; the same effect was 
found by R97a in their sample. Note that some cluster pairs 
with metallicities that differ by ~0.1 dex (Melotte 66 and 
Be 20; NGC 2141 and Be 39) have been plotted using a 
single point type and fiducial line for clarity. Clusters of 
such similar abundance are at the limit of our metallicity 
resolution. 

3 CALIBRATING THE METALLICITY SCALE 

Most previous work has found a linear relationship 
between the reduced e quiva le nt width and [ Fe/Hl 
jArmandroff fe Da Costal 119911 lOlszewski et al] I l99l| : 
Suntzeff et alJ Il99sf iDa Costa fe Armandrqfj Il9951: 
Geisler et alT Il995l iRutledee et all Il997bl: ICole et all 
20001 iTolstov et all 1200 it ISmecker-Hane et alJ 120031 
Kraft^^rvansn2003T) . The recent standard has been to 
calibrate to the relation between W' and [Fe/H] on the 
CG97 scale, as given by R97b. Given that we have derived 
W' on a totally different scale from R97b, we will need to 
rederive this relation. 

Figure 21 shows the cluster reduced equivalent widths 
plotted against the reference abundance on the CG97/F02 
scale. The best-fitting line through the data is shown, with 
the residuals to the fit shown in the bottom panel. The equa- 
tion of the fit is 

[Fe/H]gg 97 = (-2.966 ± 0.032) + (0.362 ± 0.014) W', (5) 

with r.m.s. scatter a = 0.07 dex. With this set of reference 
abundances, there is no obvious age effect on the calibration, 
nor is there any dramatic deviation from the straight-line fit 
for metallicities beyond the range of the R97b calibration. 
The only cluster that deviates by as much as la from the fit 
is Be 39; we will discuss the effects of alternate metallicity 
scales in section |S] below. 

Some curvature may be present, based on a pre- 
dicted loss of CaT index sensiti vity at high metallicity 
JPi'az. Terlevich fc Terlevichlll98Sft. Evidence fo r just such 
an effect was presented bv lcarretta et alJ i200ll) . who were 
driven primarily by upward revisions of the abundance of 
the metal-rich bulge globulars NGC 6528 and NGC 6553. 
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If we follow ICarretta et alJ (1200 if) in adding a W' 2 term to 
equation we find the coefficient of the quadratic term to 
be formally insignificant, and the scatter is not measurably 
reduced. We conclude that a linear relation suffices to reli- 
ably predict [Fe/H] on the currently adopted scale from our 
measurements of W'. However, there are two caveats: 

(i) Homogeneous high-resolution abundances for the open 
clusters are not yet available, so upward revision of their 
reference metallicities could change this conclusion. 

(ii) Age-metallicity effects could be conspiring to produce 
the observed linear relation, i.e., perhaps surface gravity 
and/or temperature effects cause a variation in W' with age 
that we do not observe because we lack the dynamic range 
in age at high metallicity to do so. 

Both of these points are further discussed below. 



4 THE OPEN CLUSTER TRUMPLER 5 

Trumpler 5 was discovered nearly 75 years go by visual in- 
specti on of pho tographic plates taken by E.E. Barnard in 
1894 jTrumr>lerll930h . The cluster lies in the Galactic plane 
» 20° from the anticentre, in a region of high and vari- 
able interstellar reddening that has complicated attempts 
to derive its fundamental parameters. Various CMD studies 
have been under taken, focusing on the red giant branch in 
the 1970's (e.g.. iKalinowski et al.lll974h . and reaching the 
unev olved main-sequence in the 1990's (K98; iKim fc Sund 
2003). Authors have converged on values near 3 kpc for the 
distance, E(B— V) = 0.6 for the reddening, and Mv = —5.8. 
K98 recognized that this makes it one of the most massive 
known open clusters in the Milky Way; we estimate it to be 
some 5 times as massive as M67, and within a factor of 10 
of the massive, intermediate-age clusters of the Magellanic 
Clouds. 

Tr 5 can thus be an important cluster for the study of 
intermediate-age stellar populations, cluster evolution, and 
population gradients in the Galaxy. Howeve r, no reliable 
metal licity measurement has yet been made. IKim fc Sund 
(2003) estimate [Fe/H] = —0.3 based on the photometry of 
K98; K98 simply took the metallicity to be equal to that of 
M67. Because of age-metallicity degeneracy, the age is simi- 
larly uncertain; the problem is worsened by the effects of red- 
dening . From the same photometry, K98 and lKim fc Suns! 
(2003) derive 4.1 Gyr and 2 .4 Gyr, respectively, wh ile the 
morphological age index of Ijanes fc Phelpsl dl994f) yields 
4.9 Gyr. 

We included Tr 5 in our observing plans for three rea- 
sons: 1) to obtain the first accurate metallicity assessment 
for this cluster; 2) to make a testable prediction on our cali- 
bration of the Ca II triplet method at high metallicities and 
intermediate ages; and 3) to lay the groundwork for subse- 
quent use of red giants in Tr 5 as metallicity standards for 
application of the CaT method from 8-10 metre class tele- 
scopes (for which nearby clusters are inconveniently bright). 

The observed stellar sample towards Tr 5 is given in 
Table^] with the adopted V magnitudes, the derived radial 
velocities, and equivalent widths. No radial velocities have 
previously been published, so we were unsure whether or not 
we could reliably discriminate between cluster members and 



the field. However, we found a well-determined mean helio- 
centric radial velocity V r = +54 ±5 km sec -1 , where the 
uncertainty reflects the r.m.s. scatter of the member stars 
about the mean. This squarely assigns Tr 5 to the Galactic 
disc, but it may be lagging in rotation by 10-20 km sec -1 
compared to the thin disc- which is unsurprising, given its 
estimated age. 

Two stars were excluded from further analysis based on 
their blueshifts of 40 and 60 km sec -1 respectively relative 
to the cluster- consistent with them being dwarfs in the so- 
lar neighborhood. We initially chose our stars based on the 
2MASS (K s , J-K s ) CMD, but later inspection of the K98 
(V, B— V) data showed that three of the stars in the south- 
ern part of the cluster appeared to be differentially reddened 
away from the main cluster RGB; these stars are identified 
and excluded from further analysis. Two of the stars near 
the observed tip of the red giant branch were treated with 
care because the cluster is known to harbor M giants, for 
which TiO veiling can bias the CaT-derived metallicities 
dOlszewski et al.lll99ll) . However, any TiO bands must be 
weak, or they would have been obvious in our spectra. We 
decided to include the 2 possible early M stars because they 
don't significantly change the cluster metallicity, so any con- 
tinuum veiling bias must be small. 

The equivalent widths and magnitudes of the ten re- 
maining RGB stars are shown in Figure^] From these mea- 
surements and Equation[2]we derive the metallicity of Trum- 
pler 5 [Fe/H] = —0.56 ±0.11. Lines corresponding to this 
value and range are shown in the figure. This is lower than 
the metallicities estimated from the K98 CMDs, but consis- 
tent with its large galactocentric radius (~ 11 kpc). A full 
re-analysis of the CMDs is b eyond the scope of th is paper, 
but we use isochrones from iGirardi et alJ ll2000l) interpo- 
lated to the correct metallicity to rederive the cluster age. 
Using the VI data, we find good agreement with the main- 
sequence, subgiant branch, and lower RGB for a distance of 
2.8 kpc, a reddening E(B— V) = 0.66 mag, and an age of 3.0 
±0. 5 Gyr. The age is b etween the values reported by K98 
and lKim fc Sund 12003). for 10-14% higher reddening and 
7-1 8% smaller distanc e. 

IKalinowski et all dl974h observed that the carbon star 
V493 Mon is projected within 2'5 of the cluster centre and 
has typical optical magnitudes for an N-type C star if as- 
sumed to lie at the cluster distance and reddening. This is 
also true if one conside rs the near-infrared photometry from 
lliiima fc Ishidal (I1978T) . Assigning a membership probabil- 
ity will require precise radial velocities and proper motions; 
however, a highly discrepant radial velocity is sufficient to 
rule out membership. We find V r = +46 km sec -1 , close to 
the cluster mean. Considering its projected distance within 
1/4 of the cluster radius of the centre, the surface density 
of field carbon stars, the optical photometry, the infrared 
photometry, and the radial velocity, V493 Mon seems to be 
a more than likely member of Trumpler 5. 

5 APPLICATION TO COMPOSITE 
POPULATIONS 

The CaT method is the best available route to metallic- 
ity measurements of RGB stars in ne arby galaxies (e.g., 
ISuntzeff etaDll993l : lT"olstov et al.|l200ll) . There are two po- 
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tential difficulties in its application. First, the reduced equiv- 
alent width, and its correlation with [Fe/H], are only defined 
and calibrated to stars older than ~10 Gyr and more metal- 
poor than ~ | solar. While this does not present severe 
problems for some of the smallest dwarf spheroidals (Draco, 
Ursa Minor, Tucana), most Local Group galaxies contain 
substantial populations that violate either the age or metal- 
licity range of the R97b calibration We have 

begun to address this calibration issue with these data. 

Second, star clusters give a strongly biased view of the 
star-formation history of a galaxy, so the field stars are of 
primary interest. Yet in a composite sample, the unknown 
mixture of ages and metallicities makes it impossible to as- 
sign a unique horizontal branch magnitude that is to be used 
in the abundance determination of an individual star. The 
observed horizontal branch (or more generally, red clump) 
is more extended in the field than the cluster, because of the 
variation of clump magnitude with age and metallicity (e.g., 
IColdll998trGirardi fc Salarisll200ll) . 

A first approach to this problem was made by 
iDa Costa fc Hatzidimitrioul il998h in their study of SMC 
clusters. These authors measured the mean red clump mag- 
nitude in each cluster to create their reduced equivalent 
widths, and then applied a theoretically-derived correction 
factor to the metallicity, based on the variation of clump 
magnitude with age. The corrections were of order —0.05 
dex for the age range from 3.5-10 Gyr, smaller than the 
precision of the abundances. We dispense with this type of 
correction here, because we have included younger clusters 
in our calibrations, and the variation of their red clump mag- 
nitu de is thus f ol ded in to Equation 

ICole et all feOOOT) estimated the bias introduced 
by the use of a composite Vhb in the calculation 
of W' for field stars, and fo und similar results to 
IDa Costa fc Hatzidimitrioul <ll998t) . We make a new, em- 
pirical measurement of the bias without recourse to the- 
oretical models by considering all of our star clusters to- 
gether ("mixed"), rather than individually ("unmixed"). We 
do this by shifting each cluster to a common apparent dis- 
tance modulus (m— M)v, and calculating the mean horizon- 
tal branch magnitude of the mixed population, (My) = 0.63, 
with an r.m.s. scatter of 0.14 mag. Using new values of 
V— Vhb calculated with this mean HB magnitude, we re- 
computed the metallicities of each cluster. 

The result of this exercise is shown in Figure |5] The dif- 
ference between the mixed vs. the unmixed metallicities is 
plotted against the unmixed metallicities for the individual 
stars in each cluster. Since the difference is constant within a 
cluster, the clusters spread out along horizontal lines in this 
plot. The errorbars in the y-coordinate give the r.m.s. spread 
in metallicity of each cluster. Figure|H]clearly shows that the 
bias introduced by variation of horizontal branch/red clump 
magnitude with age is smaller than the intrinsic precision of 
the method. At low abundances, a trend with metallicity 
is apparent, but at the high end of the range, the effect is 
not clear, perhaps due to errors in the distances of the clus- 
ters, or perhaps due to a complicated interaction between 
age and metallicity effects. Figure [U is empirically support- 
ive of the theoretically-mot i vated corre ctions sugge s ted b y 
IDa Costa fc Hatzidimitrioul <1998ft and ICole et all fcOOd) . 
The small size of the offsets indicates that application of the 
CaT method to field populations of mixed age and metallic- 



ity does not introduce large biases in the derived metallici- 
ties, within the range of parameters considered here. 



6 SUMMARY & DISCUSSION 

In connexion with a programme to determine the age- 
met allicity relation a nd radial abundance gradient of field 
red g iants in the LMC llCole et all2000tlSmecker-Hane etafl 
2003, Cole et al. 2003 in preparation), we have obtained 
medium-resolution spectra of red giants in 11 Galactic 
star clusters at the near-infrared Ca II triplet (CaT) us- 
ing the FORS2 multiobject spectrograph at the 8.2-metre 
VLT/UT4. Our spectra have extremely high signal-to- noise, 
ranging from 20 < S/N < 85. Following many previous au- 
thors (e.g., R97b) we have derived a linear relation between 
the equivalent width of the CaT and the cluster metallicities. 

For the first time, our sample includes open clus- 
ters with a range of metallicities from —0.6 < [Fe/H] 
< -0.2 and ages from 2.5 < (age/Gyr) < 7~This al- 
lows an initial assessment of the effects of age on the 
metallicity measurement via the CaT. While there is is 
substantial scatter and the sample size is small, we find 
no evidence for a decrease in sensitivity of the CaT at 
young ages and high metallicities. This behavior is theoreti- 
cally sup ported by non-LTE stellar atmosphere calculations 
(e.g.. IjQrgensen et aljll992T) and is in good agreement with 
comprehensive empirical calibrations of the CaT behavior 
over a wide range o f stellar atmospheric parameters (e.g., 
ICenarro et all |2002|) . However, our results disagree with 
earlier results based on more limited stellar samples (e.g., 
|Joi^s u Allom_fc_Joj^^ fc Picklesl 

Il986t iDfaz. Terlevich fc Terlevichlll989l) . The fact that we 
see no loss of metallicity sensitivity above [Fe/H] « —0.5 is 
interesting in light of the observed metallicity independence 
of calcium line strength in the integrated light of clu sters and 
galaxies above this level fe.g.. IVazdekis et al.ll2003T) . 

One difference between our analysis and earlier work 
is that we fit a composite function- the sum of a Gaussian 
and a Lorentzian- to the line profiles. This is mandated by 
the line shapes, which indicate that for abundance [Fe/H] 
> —0.6, the damping wings of the lines are so strong and 
broad that we begin to resolve the line even at our resolution 
of ~2.5 A. 

Using our data, we obtain a linear relation between 
the reduced equivalent width of a cluster, W', and metal- 
licity on a scale formed by the high-dispersion results of 
CG97 for the globular clusters, and medium-dispersion Fe 
indices from F02 for the open clusters. There is no indi- 
cation that a second -order term in W', as advocated by 
ICarretta et al J i200lfl . is needed to fit the data. However, our 
sample does not include old, high-metallicity clusters (such 
as bulge globulars NGC 6528 and NGC 6553), clusters of su- 
persolar metallicity (such as NGC 6791 and NGC 6253), or 
young clusters of very low abundance (such as the clusters 
of the Small Magellanic Cloud) . Until such time as the age- 
metallicity plane is filled, we cannot definitively rule out an 
age effect on the calibration. Filling in the high-metallicity 
regime at a large range of ages will have important con- 
sequences for the determination of abundances in the M31 
halo. 

An important point is that the globular cluster metallic- 
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ity scale is not absolutely established. The CG97 scale differs 
in a nonlinear way from the average of spectrophotometric 
indices compiled bv lZinn fc Westl (^984) . In turn, CG97 may 
be superseded by abundance scales derived from Fe II lines, 
rather than Fe I lines, beca use of the newly-reco gnized im- 
portance of non-LTE effects ijKraft fc Ivansl2 bo3). The open 
cluster abundances are even more uncertain, and will remain 
so until large-scale, homogeneous surveys for high-dispersion 
spectra can be undertaken. This type of analysis is presently 
available for only a small number of stars in a small number 
of clusters; among our sample, Melotte 66 and M67 have 
been analyzed. Photometric metallicity indic ators based on 
DDO -system or Stromgren photometry (e.g.. iTwarog et"al] 
Il997l) often give quite different results from the indices tab- 
ulated by Friel and collaborators. The choice of metallicity 
scales will of course have a strong effect on the calibration 
of the CaT. We tabulate some alternate cluster metallicities 
in Table El Columns 1-4 give the cluster name, W derived 
from Equation^ the reference abundance from Tabled and 
the derived metallicity using Equation^] Alternate metallic- 
ities are given in column 5, sorted by origin, and referenced 
in the last column. Using the tabulated values of W and 
[Fe/H] a i t , we derive the following relation (excluding Tr 5): 

[Fe/H] a it = -2.629 + 0.058 W' + 0.037 W' 2 . (6) 

Significant curvature is introduced into the relationship be- 
cause the globular cluster metallicities have been shifted 
down, while the open cluster metallicities, especially 
NGC 2141 and M67, have been shifted upwards. The clus- 
ter metallicities derived from Equation HJ are given in Table 
Hj] We cannot comment on the correct choice of metallicity 
scales, but we note that although the linearity of the CaT 
calibration seems to be a fortuitous result of one such choice, 
the one-to-one correspondance between W' and [Fc/H] holds 
no matt er the choice of refe rence metallicities (as already 
noted bv lKraft fc Ivansll2003l) . 

We have measured the bias introduced by taking the 
mean horizontal branch/red clump magnitude of a compos- 
ite field population in determining the quantity V— Vhb- 
Shifting all 10 clusters to a common apparent distance mod- 
ulus and taking the mean Vhb of the mixed population re- 
sults in derived metallicities that shift by ±0.05 dex relative 
to the unmixed results. This is negligible within the uncer- 
tainties, and shows that measurements of individual field 
red giants yield reliable metallicities despite the lack of a 
unique Vhb in that scenario- although the additional uncer- 
tainty should be considered in the error analyses. The mixed- 
population bias is clearly not random at the low-metallicity 
end of our calibration sample, where the age range is small. 
This is consistent with theoretical expectations for the var i- 
ation of Vhb with metallicity fe.g.. lGirardi fc SalarisllioOll) . 
At the high-metallicity end, the picture becomes more com- 
plicated, with age and abundance effects combining with 
distance uncertainties to introduce additional scatter. 

We have made the first measurement of the metal- 
licity of the large, old open cluster Trumpler 5, finding 
[Fe/H] = —0.56 ±0.11. While interesting in itself, this also 
represents a testable prediction of our calibration of the 
CaT. With our metallicity estimate, the most likely age 
for Tr 5 is 3.0 ±0.5 Gyr. The relatively low metallicity of 
Tr 5 places the cluster in a category with objects such as 
Melotte 66. It somewhat diminishes the utility of Tr 5 as an 



intermediate-age, high-metallicity calibrator for the Ca II 
triplet, but makes it easier to theoretically understand the 
possible presence of a carbon star in the cluster. The carbon 
star, V493 Mon, has a radial velocity within 8 km sec -1 of 
the cluster mean. Given the large number of independent 
data consistent with membership, it would be surprising if 
V493 Mon turned out to be a field interloper. However, more 
precise velocities, as well as proper motions, are needed to 
assign V493 Mon to Tr 5 for certain. This raises the excit- 
ing possibility that detailed study of the cluster will yield its 
turnoff mass (« 1.141q 04 Mq), and therefore a direct mea- 
surement of the mass of the progenitor of an intermediate- 
mass, moderate-metallicity carbon star. 

We are further refining the empirical calibration of 
the CaT by enlarging the sample of high-metallicity, 
intermediate-age cluster giants with high signal-to-noise 
spectra. TLB and TSH have measured the spectra of nu- 
merous giants in several clusters, at both high- and low- 
dispersion, using the Shane 3-metre Telescope at Lick Obser- 
vatory. The cluster sample spans the abundance range from 
-2.3 < [Fe/H] < ±0.4 and ages > 2 Gyr. In addition to the 
new CaT measurements, we will determine [Ca/H] abun- 
dances for the clusters using model atmosphere analyses of 
the strengths of the plentiful neutral calcium lines (which are 
theoretically much better understood) in the high-dispersion 
spectra. We will calibrate the relationship between the CaT 
index and [Ca/H], eliminating the need to assume a specific 
[Ca/Fe] vs. [Fe/H] relationship. This is an important step 
for application of the CaT technique across diverse environ- 
ments which may have experienced very diffe rent histories o f 
chemical evolution (e.g., dwarf galaxies iTolstov et alT2 003y 

Our results are a striking reaffirmation of the util- 
ity of the CaT for abundance measurements of red giants 
in intermediate-age field populations. Red giants are the 
brightest common stars among old stellar populations, and 
so they are currently the only practical targets for abun- 
dance analyses of stars aged as 2-14 Gyr in Local Group 
galaxies beyond the Milky Way halo. In these systems, the 
typical magnitude of the RGB tip is I > 20, and measure- 
ment of the CaT from low-resolution spectra is the most 
direct available method to accrue the large samples of abun- 
dance measurements needed to assess the chemical evolution 
of systems beyond the Milky Way halo. 
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Table 1. The Cluster Sample 



Cluster 


a (J2000) 


S (J2000) 


[Fe/H] 


Age (Gyr) 


(m-M)v 


V HB 


Reference 


NGC 104 (47 Tuc) 


00:26:33 


-71:51 


-0.70 ±0.07 


10.7 


13.37 


14.06 


1,2,3 


NGC 1851 


05:14:15 


-40:04 


-0.98 ±0.06 


9.2 


15.47 


16.09 


2,3,4 


NGC 1904 (M79) 


05:24:12 


-24:31 


-1.37 ±0.01 


11.7 


15.59 


16.15 


1,2,3 


NGC 2298 


06:48:59 


-36:00 


-1.74 ±0.06 


12.6 


15.59 


16.11 


1,2,3 


NGC 4590 (M68) 


12:39:28 


-26:45 


-1.99 ±0.10 


11.2 


15.19 


15.68 


1,2,3 


Berkeley 20 


05:32:34 


+00:10 


-0.61 ±0.14 


4.9 


15.01 


15.70 


5,6 


NGC 2141 


06:03:00 


+10:30 


-0.33 ±0.10 


2.8 


14.16 


14.90 


5,7 


Melotte 66 


07:26:28 


-47:41 


-0.47 ±0.09 


6.3 


13.63 


14.52 


5,8 


Berkeley 39 


07:46:45 


-04:41 


-0.26 ±0.09 


7.2 


13.78 


14.28 


5,8 


NGC 2682 (M67) 


08:51:25 


+11:48 


-0.15 ±0.05 


6.3 


9.98 


10.54 


5,8 


Trumpler 5 


06:36:34 


+09:27 


-0.3: 


4 : 


14.40 


15.13 


9,10 



Refere nces to metallicity, distance mod ulus, and Vhr va lues: (l_) JC^£rejAa_j£Grattonl ll997[): (2) Salaris fc Weissj l|2002[) : (31 [ Harris! 
<1996ft : (4l lRutledge et alj Jl997bl) : C5) lFriel et alTll2002l); ( 6) lMacMinn et alj 1994TT (7) iRosviclj <1995l) : (8) ISaraiediniHl999l) : (91 

iKaluznvl lll998h : Clo HKim fc Sunel <200fl . 
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Table 2. Observed Stars: Globular Clusters 



Star 


V 


V r (km s" 1 ) 


SW (A) 


47 Tuc: Lee (19771 










L5309 


12.21 


-21.8 


±7.6 


7.98 ±0.10 


L5310 1 


13.41 


9.1 ±7.6 


6.46 ±0.10 


L5312 


12.18 


-12.3 ±7.5 


7.91 ±0.10 


L5418 


15.33 


-20.6 ±7.9 


5.77 ±0.16 


L5419 


14.05 


-22.2 ±7.8 


6.13 ±0.12 


L5422 


12.47 


-23.1 ±7.6 


7.38 ±0.10 


L5527 


13.56 


-25.5 ±7.6 


6.79 ±0.09 


L5528 


14.40 


-28.5 ±7.6 


6.37 ±0.18 


L5530 


13.14 


-18.6 ±7.6 


6.92 ±0.08 


' Radial velocity nonmember, excluded from analysis. 


NGC 


1851: Stetson (19811 








003 


13.60 


324.5 


±7.4 


7.52 


±0. 


.15 


065 


15.75 


322.1 


±7.5 


5.67 


±0. 


.08 


095 


13.57 


334.1 


±7.4 


6.94 


±0. 


.11 


107 


14.50 


330.7 


±7.5 


6.78 


±0 


.13 


109 


14.85 


334.3 


±7.4 


5.96 


±0. 


.10 


112 


13.80 


327.3 


±7.4 


6.85 


±0 


.11 


123 


16.21 


329.7 


±7.4 


4.99 


±0. 


.08 


126 


14.29 


321.5 


±7.5 


6.95 


±0. 


.11 


160 


15.51 


326.4 


±7.7 


5.85 


±0. 


.09 


175 


16.22 


332.3 


±7.9 


6.07 


±0. 


.10 


179 


16.45 


320.4 


±7.8 


5.66 


±0. 


.08 


195 


15.82 


327.6 


±7.5 


5.75 


±0. 


.09 


209 


14.08 


333.4 


±7.6 


7.05 


±0. 


.20 


231 


15.93 


332.8 


±7.8 


5.89 


±0. 


.10 


275 


14.98 


333.8 


±7.6 


6.79 


±0. 


.13 


NGC 


1904: Stetson & Harris (19771 








006 


15.27 


203.2 


±7.1 


4.91 


±0. 


.08 


011 


15.74 


213.5 


±7.5 


4.85 


±0. 


.10 


015 


13.22 


204.8 


±7.5 


6.67 


±0. 


.10 


045 


15.58 


206.8 


±7.8 


4.95 


±0. 


.08 


089 


14.71 


206.5 


±7.8 


5.51 


±0. 


.09 


091 


16.25 


207.6 


±7.8 


4.76 


±0. 


.08 


111 


15.62 


206.9 


±7.6 


5.18 


±0. 


.08 


115 


15.95 


211.0 


±7.7 


4.65 


±0. 


.10 


138 


16.16 


206.3 


±7.8 


4.36 


±0. 


.08 


153 


13.44 


206.1 


±7.4 


6.47 


±0. 


.11 


160 


13.02 


200.3 


±7.5 


6.42 


±0. 


.10 


161 


15.86 


216.5 


±7.8 


4.05 


±0. 


.07 


176 


14.95 


215.9 


±7.6 


4.74 


±0. 


.08 


209 


15.02 


210.8 


±7.5 


5.63 


±0. 


.11 


224 


16.09 


220.4 


±7.9 


4.63 


±0. 


.09 


237 


14.19 


220.8 


±7.5 


5.77 


±0. 


.11 


241 


13.61 


228.5 


±7.5 


5.80 


±0. 


.09 



Table 2 - continued Observed Stars: Globular Clusters 



Star 


V 


V r (km s" 1 ) 


SW (A) 


NGC 


2298: Alcaino & Liller (19861 




6 


13.82 


156.5 ±7.6 


4.97 ±0.07 


12 


14.23 


146.6 ±7.5 


4.84 ±0.07 


15 


14.86 


151.2 ±7.6 


4.31 ±0.06 


22 


15.30 


162.0 ±7.4 


4.06 ±0.06 


25 


15.69 


153.2 ±7.5 


4.15 ±0.06 


S156* 


15.53 


148.5 ±7.6 


4.14 ±0.05 


S172* 


15.61 


160.6 ±7.7 


3.97 ±0.05 


* ID and magnitude from T. Smecker-Hane, unpublished. 


NGC 4590: Harris (19751 




1-2 


14.95 


-96.0 ±7.6 


2.98 ±0.04 


1-49 


14.62 


-93.6 ±7.8 


3.04 ±0.05 


1-74 


14.59 


-103.7 ±7.7 


3.21 ±0.05 


1-119 


13.62 


-89.1 ±7.6 


3.78 ±0.06 


1-239 


14.19 


-86.0 ±7.6 


3.63 ±0.05 


1-256 


12.64 


-90.9 ±7.6 


5.20 ±0.08 


1-258 


14.24 


-90.2 ±7.7 


3.14 ±0.05 


1-260 


12.52 


-91.4 ±7.5 


4.85 ±0.07 


11-47 


15.03 


-85.8 ±7.8 


3.19 ±0.05 
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Table 3. Observed Stars: Open Clusters 



Star V V r (km s" 1 ) SW (A) 



Berkeley 20: MacMinn et al. ('1994) 



05 


14.80 


73.9 ±7.5 


7.56 ±0.12 


08 


15.15 


79.2 ±7.5 


7.41 ±0.12 


12 


16.21 


83.8 ±7.5 


6.61 ±0.12 


22 


16.90 


74.6 ±7.5 


6.19 ±0.13 


29* 


17.14 


42.1 ±7.7 


5.52 ±0.17 



t Radial velocity nonmember, excluded from analysis. 

NGC 2141: Burkhead. Burgess fc Haisch (1972) : Rosvick (19951 



1-3-21 


14. 


.18 


31.8 ±7.5 


8.04 ±0.18 


1-4-05 


14. 


.61 


50.7 ±7.5 


7.58 ±0.13 


3-2-18 


13. 


.05 


26.0 ±7.5 


9.35 ±0.18 


3-2-34 


14. 


.95 


27.6 ±7.5 


7.26 ±0.17 


3-2-40 


13. 


.25 


28.9 ±7.4 


8.80 ±0.17 


3-2-52 


14. 


.36 


33.2 ±7.5 


7.67 ±0.14 


3-2-56* 


15. 


.54 


-11.3 ±7.5 


4.15 ±0.09 


4-08 


14. 


.80 


23.2 ±7.6 


7.76 ±0.13 


4-09 


13. 


.27 


28.7 ±7.6 


8.78 ±0.14 


4-13 


15 


.15 


32.1 ±7.4 


6.79 ±0.14 


4-14 


15. 


.52 


28.8 ±7.4 


7.38 ±0.12 


4-24 


14. 


.77 


31.4 ±7.5 


7.64 ±0.13 


4-25 


14. 


.13 


31.3 ±7.5 


8.11 ±0.18 


5-09 


14. 


.62 


33.0 ±7.5 


7.69 ±0.14 


5-13 


13. 


.96 


32.0 ±7.5 


7.92 ±0.15 


5-13a 


15. 


.00 


49.3 ±7.4 


7.00 ±0.16 


* Radial velocity nonmember, 


excluded from analysis. 



Melotte 66: Anthonv-Twaroe. Twaroe * 


r. McClure (19791 


1205 


14.18 


18.3 ±7.5 


6.86 ±0.11 


2107 


14.70 


14.9 ±7.5 


6.57 ±0.11 


2133 


13.20 


10.9 ±7.5 


7.91 ±0.12 


2226 


14.08 


18.9 ±7.5 


7.11 ±0.19 


2233 


15.44 


14.4 ±7.4 


6.51 ±0.09 


2244 


14.47 


16.4 ±7.6 


6.97 ±0.09 


2261 


13.64 


41.2 ±7.5 


7.46 ±0.16 


3101 


14.98 


15.6 ±7.7 


5.03 ±0.10 


3133 


14.08 


17.2 ±7.6 


6.98 ±0.12 


3229 


13.93 


7.7 ±7.5 


7.08 ±0.31 


3235 


14.66 


18.3 ±7.7 


6.89 ±0.10 


3260 


14.56 


16.2 ±7.4 


6.90 ±0.12 


4151 


12.69 


12.7 ±7.6 


8.23 ±0.13 


4265 


14.07 


7.4 ±7.5 


6.97 ±0.12 



Table 3 - continued Observed Stars: Open Clusters 



Star V Vrfkmr 1 ) SW (A) 



Berkeley 39: Kalu znv fc Richtler (1989) 



UUz 


1 q no 

to.uy 


00.0 ± i .O 


Q QO _i_n 1 c 


UUo 


lo.OO 


01. U ± 1 .4 


Q Q£ _i_n 1 Q 


UUO 


lo.oo 


00.0 ± i .0 


i fiQ _i_n i a 
(.Do ztU.lu 


nno 
uuy 


1 A 0*3 


£fl "7 -1-7 K 

oo./ ±t.o 


(.14 ztU.lo 


012 


14.34 


55.9 ±7.4 


7.37 ±U.14 


013 


14.35 


58.2 ±7.5 


6.87 ±0.16 


016 


14.48 


61.4 ±7.5 


7.28 ±0.15 


017 


14.75 


54.0 ±7.6 


6.72 ±0.12 


018 


14.76 


57.9 ±7.5 


6.88 ±0.14 


028 


15.40 


58.9 ±7.5 


6.60 ±0.15 


M67*: 


Johnson <! 


^Sandage (19551: Sanders 


(19771 


F104 


11.20 


33.5 


7.13 ±0.16 


F105 


10.30 


34.3 


8.00 ±0.16 


F108 


9.72 


34.7 


8.36 ±0.17 


F135 


11.44 


34.3 


7.10 ±0.13 


F141 


10.48 


33.6 


7.73 ±0.14 


F164 


10.55 


33.3 


7.40 ±0.15 


F170 


9.69 


34.3 


8.28 ±0.19 


1264§ 


11.74 


34.9 


5.83 ±0.23 



* Radial velocities from i Matliievi^^al i (1986). 

§ Blend (Monteomcrv, Marschall & Janes 1993); excluded from analysis. 
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Table 4. Observed Stars: Trumpler 5 


Start 


V (mag) 


V r (km s" 1 ) 


SW (A) 


Note 


JaO / UO 


14 54 


-o.y in / .0 


7 71 _|_n ITi 




JaOOO 1 * 


14 48 


r:o 7 _l_7 tr: 






-LVOUUU 


1 A 


^fl 7 4-7 A 


7 04 -1-0 




K2579 


15.59 


47.8 ±7.5 


6.37 ±0.03 




K2565 


15.79 


60.5 ±7.5 


6.12 ±0.03 




K2324 


15.03 


53.1 ±7.5 


6.70 ±0.02 




K2280 


15.16 


47.8 ±7.4 


6.54 ±0.03 




K1935 


12.87 


53.5 ±7.5 


8.63 ±0.03 


2 


K1834 


16.24 


50.8 ±7.7 


6.21 ±0.05 




K1401 


15.29 


14.8 ±7.5 


6.60 ±0.03 


1 


K1305 


14.24 


64.9 ±7.5 


8.12 ±0.03 


3 


K1277 


14.92 


46.3 ±7.2 




1 


K1214 


13.91 


55.7 ±7.5 


7.61 ±0.03 


2 


K1026 


14.73 


52.0 ±7.3 


6.48 ±0.10 




K833 


14.91 


54.4 ±7.4 


6.53 ±0.04 


3 


K488 


16.51 


57.1 ±7.7 


7.31 ±0.05 


3 



t From Table 3 of Kaluzny (1998). 

Notes: (1) radial velocity nonmember; (2) weak TiO?; (3) Prob- 
ably differentially reddened; (4) Carbon star = V493 Mon. 



Table 5. Cluster Radial Velocities 







<r(V r ) 


prev. V r 




Cluster 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


Reference 


47 Tuc 


-22 


5 


-18.7 


1 


NGC 1851 


329 


5 


320.5 


1 


NGC 1904 


211 


7 


206.0 


1 


NGC 2298 


154 


6 


148.9 


1 


NGC 4590 


-92 


6 


-94.3 


1 


Be 20 


78 


5 


70 


2 


NGC 2141 


33 


5 


64 


3 


Mel 66 


16 


8 


23 


1 


Be 39 


58 


2 


55 


2 


M67 






33 


2 


Tr 5 


54 


5 







References: (11 iHarrisl l ll996ft: (21 IFriel et alJ fc002l) : (3) 
iFriel. Liu, Janea <19S9fr : m TFriel fc JanesTTr993T . 
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Table 6. Cluster Metallicity Data 



W [Fe/H] [Fe/H] [Fe/H]t [Fe/H] 

Cluster (A) (ref) (Eqn.|SJ (alt.) (Eqn.|HJ Reference 



47 Tuc 


6.36 


±0 


03 


-0.70 


±0.07 


-0 


66 


±0 


09 


-0.88, -0.71 


-0 


77 


±0 


13 


1,2 


NGC 1851 


5.55 


±0 


02 


-0.98 


±0.06 


-0 


96 


±0 


12 


-1.19, -1.33 


-1 


17 


±0 


16 


1,2 


NGC 1904 


4.40 


±0 


02 


-1.37 


±0.01 


-1 


37 


±0 


11 


-1.64, -1.68 


-1 


66 


±0 


12 


1,2 


NGC 2298 


3.54 


±0 


03 


-1.74 


±0.06 


-1 


69 


±0 


07 


-2.07, —1.81 


-1 


96 


±0 


06 


1,2 


NGC 4590 


2.48 


±0 


03 


-1.99 


±0.10 


-2 


07 


±0 


09 


-2.43, -2.09 


-2 


25 


±0 


26 


1,2 


Be 20 


6.91 


±0 


03 


-0.61 


±0.14 


-0 


47 


±0 


07 


-0.4 


-0 


17 


±0 


11 


3 


NGC 2141 


7.47 
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-0 
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±0 
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±0 


03 
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-0.05, 0.00 
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±0 


09 
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±0 


03 






-0 


56 


±0 


11 


-0.3 


-0 


62 


±0 
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7 



' Values in roman type are from high-disper sion spectroscopy; v alues in it alic type are from phot omet ric or spectrop hotometric indices. 
References: (1) Kraft fc Ivans! i2003ft : f2~)IZinn Wesj <1984h: (?,) iDurgapal et alJ fcOOlft: (4) iTwarog et all Il997tl : (5) 
Gratton fc Contarinil Il994l : (61 IShetrone fc SandauisJ l200Cll : f7') lKim fc Sund l200a> ! 
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Figure 1. Sample spectra for stars in three clusters, representative of the range of line strengths in our targets. The spectra have been 
sky-subtracted, continuum-normalized, and Doppler-shifted to the rest frame for comparison. Each spectrum is offset by a constant 
amount from the one below. Stars are identified by cluster name and metallicity from Table and ID from Tables [5] and [3] Most of the 
"noise" in the most metal-rich spectrum is due to the increasing strength of a large number of relatively weak atomic lines, mainly Fe I. 
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Figure 2. (Top panel) Summed equivalent widths of the 3 CaT lines measured using a Gaussian fit, compared to the same lines fit with 
the sum of a Gaussian and a Lorentzian function. The dashed line shows the best linear fit to the globular cluster stars. (Bottom panel) 
The residuals to the fit are plotted, showing how a Gaussian fit increasingly fails to measure the rising fraction of absorption in the far 
line wings as the lines get stronger. 
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Figure 3. The CaT equivalent widths of RGB stars in 10 star clusters compared to their magnitude above their respective cluster 
horizontal branches. The solid lines show the best-fitting lines to each cluster, assuming the best common slope of 0.73 A mag . Clusters 
are represented by: closed squares (NGC 4590), open squares (NGC 2298), closed triangles (NGC 1904), open triangles (NGC 1851), 
crosses (47 Tuc), stars (Mel 66 and Be 20), open circles (NGC 2141 and Be 39), and closed circles (M67). The fit to NGC 4590 adopting 
the slope from R97a of 0.64 A mag -1 is shown for comparison (dotted line). Errorbars are omitted for clarity, but the typical la errors 
are shown at lower left. 
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Figure 4. (Top panel) The reference metallicity of the 10 clusters in Table HI plotted against their reduced CaT equivalent width. The 
best linear fit to the data is shown. (Bottom panel) the residuals to the linear fit are plotted; the errorbars are la errors on the reference 
abundance. 
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Figure 5. CaT equivalent widths vs. V magnitude for 10 RGB stars in the old open cluster Trumpler 5. The lines corresponding to 
the best-fitting mctallicity according to our calibration (solid), and the formal lcr rms scatter (dashed) arc shown, with reference lines 
(dotted) corresponding to —0.2 and —1.0 dex. 
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Figure 6. The difference between "mixed" and "unmixed" abundance estimates for the 10 well-studied star clusters. The "mixed" 
estimates are made by shifting all clusters to a common (m— M.)y and using the combined average red clump + horizontal branch 
magnitude to define Vhb- The errorbars show the rms scatter of the "unmixed" abundances of each cluster. 



